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ABSTRACT 

TlK overall goal of the "Heavy Oas Project- i> to develop an improved computer 
model for tlie ptedicrion of iieavier-llian-air gas dispersion in complex leirain. 
Three basic areas of investigation comprise the Project: 1) experimental - zero 
wind releases; 2) experimental - releases into wind; and 3) mathematical 
modelUng. In this paper we present our wortc on Part 1. More specifically, the 
spreading rate of a heavier-than-air gas cloud into a zero-wind environment is 
studied using now visualizaiion techniques. Smoke was used to follow the 
instantaneous release of a heavy gas onto firsUy, a smooth level surface, and 
subsequenUy onto surfaces with obstacles arranged in different configurations. 
Overhead and side-view cameras were used tp photograph the spill events, 
sun photos were obtained from fteeze-frame analysis of die video, and were 
then used to produce data on cloud fhmt distance versus elapsed time. 
Photos of the release showed a very uniform circular cloud spread onto a 
smoodi surface. A release into obstacles appeared to start out initially as a 
uniform circular cloud, but at later times die cloud developed into a diamond or 
hexagonal shape, depending on the obstacle arrangement. Surprisingly, the 
cloud spreading rate was higher widi obstacles present than for spreading on a 
smooth surface. This, however, was true only for cloud spreading through an 
open coiridor in a uniform array of obstacles. A modified MOE model 
predicted die cloud spreading history with reasonable results. 

Keywords: heavy gas. still-air. flow visualization. smooUi. obstacles, model 
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I. INTRODUCTION 

The widcspFcad Iraiupon and storage of loxic and/or inflaitinuble gases 
has increased interest in recent years in the behaviour of heavier-than-air gases 
released into the ambient atmosphere. In the event of a spill of a toxic or 
flammable heavy gas, an emergency response team must be able to respond 
quickly and efficiently to ensure that danger to the local community ii 
minimized. This requires a real-time ability to predict the flow path and 
concentration paiiem of the gas cloud for the duration of the spill event. In 
many cases, the gas is heavier than air, especially during the early stages of 
the release. The behaviour of the heavier-lhan air gas cloud depends on 
several variables such as local terrain roughness and slope, gas density, physical 
nature of the gas release (instantaneous, semi-continuous or continuous), spill 
rate and/or spill volume as well as weather conditions including wind speed and 
direction, atmospheric stability, precipitation, and humidity. 

Perhaps Uie most important parameters In determining the behaviour of the 
gas cloud are gas volume and density, wind speed and direction and local 
terrain characteristics. The fact that weather conditions play a large role in 
heavy gas behaviour enqihasizes die complex nature of die problem. 
Nevertheless, a large body of research has been, and is being, completed in an 
effort to help provide a quicic and efficient response to a spill of toxic or 
inflammable heavy gases. Much of this work has been sponsored by indusDiy 
and government in a approach similar to the joint venture at Thomey IsUnd 
(McQuaid, 198S). 

Research into heavier-Uian-air gas dispersion has proceeded on three 
fronts: 1) full-scale field trials, 2) wind tunnel simulations, and 3) theoretical 
prediction or computer simulation. Full-scale tests represent an enormous effort 
in manpower and funding. Wind tiuinel and computer simulations are 
considerably less expensive and are very useful for validation purposes. In diis 
study, the inslanuuieous release of a volume of heavy gas is simulated under 
zero-wind conditions and several different obstacle configurations. Similar tests 
are cunenlly being carried out with non-zero wind, and these will be Ute 
subject of a later paper. Ultimately, the results of flow visualization and 



111 






■ ■•^5 



'3^.- 



conccniralion moniioring from ihese lesu will be used lo improve mathemalical 
models of heavy gas dispersion. 

In the first series of lieavy gas dispersion experimenis, a Icnown volume of 
gas was released inlo a calm environment under a variety of obstacle arrays. 
For eacli of tliese releases, ihe flow of the gas cloud was observed and 
recorded on video tape. By varying Ihe size and density of the obstacles as 
well as the initial cloud dimensions, the effects of these obstacles could be 
studied and then compared to simple box model predictions of the horizontal 
cloud dimensions. Improved mathematical expressions are proposed to describe 
Ihe horizontal nxition of Ihe gas cloud in its early stages of dispersion. 

2. EXPERIMENTAL PROCEDURES 

The scale modelling of a heavy gas release consists of holding certain 
physical parameters constant for both the full scale and laboratory tests. 
Similarity criteria include Ihe Reynolds (Re) and Froude (Fr) numbers based on 
wind velocity, Richardson number (Ri) based on wind velocity and density 
difference between gas spill and air, initial density ratio, volume and rate of 
release, wind velocity profile, surface roughness, turbulence (length scale and 
intensity) and weather stability class. The Reynolds number is given by: 



Re = LU8 /H 



(I) 



where L is a length scale (arbitrarily defined), U is local wind velocity, S, is 
the air density, and \i is the viscosity. The Richardson number is: 



Ri = gA5/83lVU^ 



m 



where g is the gravitational constant, AS == 6 - &^, and S is Ihe gas 
density. 

Not all criteria can be satisfied in practice. The Reynolds criteria cannot 

be satisfied for models at 1:300 scale because the velocity over the model 
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becomes very low, and it is difficult to obtain a good simulaiion of Ihe 
boundary layer at low speeds in die wind tunnel. Fortunately, a reasonably 
complete simulation can be obtained, even in violauon of Reynolds similitude 
(Meroney, 1982), (Meroney and Neff, 1986). This "non-exact" simulation can be 
esublished by ensuring that Ihe Richardson number, Ri, is equal al both model 
and full scale. The density ratio may be increased in the model case so diat 
Ihe velocity on Ihe model may be increased. This is done because a higher 
velocity simulation is much easier to perform. 

The velodly profile over Ihe model should be similar to the fuU scale 
wind profile. The absolute values of the velocitiea ate obtained from Ri 
similitude. The surface roughness is scaled simply by geometric scaling. The 
luibulence intensity of die wind in die along-wind direction should be similar al 
model and full scale. In partial boundary layer simuladon, some of dieac 
criteria may not be satisfied exacdy. The velocity profile and surface 
roughness criteria are relauvely easy lo saOsfy, but die length scale and 
turt)ulence intensity are more difficult However, because die turbulence 
generated during die instantaneous release of a heavy gas spill is generally 
higher dian atmospheric turbulence, exact similitude is not required. Thus, 
aldiough die simulated boundary layer lacks some turtiulent elements, die model 
spill will remain a reasonable represenudon of die full scale event 

The release mechanism developed for diis study is illustrated in Figure 1. 
Design criteria included die following: 

i) scale release volume; 

ii) geometric shape; 

iii) minimum interference from release; 

iv) leak proof; and 

v) flexible pull down mechanism. 

Two release canisters were constnicied to repieseni significantly different 
volumes of gas release. Al 1:1(X) scale, die small canister (2(KX) cm^) 
represented a Thorncy Island release of 2000 m'. The large canister (14000 
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cm') represented a 1:50 scale release. Figure 2 shows the release from the 
small canister into a zero-wind environment with no blocks. A small puff of 
smoke from the lop of the release canister signals the beginning of the 
release, as shown in Figure 2- A. Less than 100 msec after the release, the 
cylinder disappean through the floor. The cloud remains cylindrical, is sheared 
uniformly along the sides, and begins to collapse before 100 msec. Figure 4-C 
and D show ihe overhead view of the spill. It can be seen that the spread of 
the cloud is uniform in all directions, indicating that the release mechanism was 
performing satisfactorily. 

The gas used in the flow visualization study was Freon-12 mixed with 
carbon dioxide and smoke generated from mineral oil. By a careful weighing of 
a known volume of the Freon-CO,-smoke mixture, the density of the heavy gas 
used in this study was found to be 3.89 kg/m . 

The test arrays used in the initial flow visualization experiments resulted 
in 18 different test configurations which are defined in Table I. Two block 
sizes were used: a small cube with 3.5 em lidei and a medium cube of 7 cm 
sides. Block coverage was either none, half or full. A schematic of full 
coverage is given in Figure 3. Half coverage has blocks only in the semicircle 
to the south of the canister location; and north of the canister there were no 
blocks (see North arrow on photos). 

In most of the experiments using the cubes, the spacing of the blocks on 
the platform was 1:1 as shown in Figure 3. The 1:1 spacing indicates that the 
distances between the blocks (as indicated by B and D on Figure 3) is equal to 
die block width (as indicated by A and B on Figure 3). In two release series, « 
staggered spacing of 1:1 was used in which alternate rows are shifted one 
spacing to the left of the row below and above it. This spacing results in a 
1:1 spacing within all rows and between the rows and a 1:3 spacing In each 
column. A 1:3 spacing indicates 3 block widths between adjacent blocks. In the 
regular 1:3 spacing, the blocks are so spaced in both the tow and column 
directions. 
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Two ouiistcr sizes were also used in the release experiments. The small 
canister has a volume of 1995 cm^ a height (indicated by G on Figure 3) of 
12.5 cm and a diameter (F on Figure 3) of 14.1 cm. The large^ canister was 
25.1 cm in diameter and 27.6 cm in height giving a volume of 13 615 cm . 

In order to better scale the How visualization records, two white lapet 
were located on model platform: one "above" ^ U» other "below" the 
canister (e.g. see Figure 4). THese tape, were located equidistant from the 
canister centre «ul were spaced 120.9 cm apart for the small canister ^ 
192 7 cm apan for the large canister. In the later releases, the upes were 
consistenUy placed 193 cm apart; white blocks were placed at the edge of the 
disk on the horizontal and veriical axel to facUiute the location of the 
canister centre (e.g. Figure 4). 

For each obstacle array. Uie instantaneous release event was recorded on 
video upe by a camera mounted overhe«l In sevend tests, a cuner. was also 
posiuoned to give a side view of the event A digital dme display «Ided to 
Ihe video tapes is seen to the upper comer of e«:h of U» photographs 
presented in this paper (e.g. Figure 4) and gives the time in hours, minutes, 
seconds and 1/30 second incremenu. In otder to analyze the history of U« 
cloud spread, the video tape was stopped at intervals and a photographic record 
„«„ie of the freeze-f«me. THese photographs were Uien used to de«.mine the 
horizontal dimensions of the cloud at several instances after release. 

3. FLOW VISUALIZATION OF STILL AIR RELEASES 

In Uiis first series of experiments, the horizontal spread of an 
instantaneous release of heavy gas under calm-wind conditions was simulated. 
Both the large and U» small gas canister, were used in the flow visual.zanon. 
Each canister offered desirable feature, in «ldition to the simuladon of a 
number of different source/obstacle ratios. The large canister produced a 
more visible cloud than die small cylinder which faciliuted tracking of the 
cloud front. The disadvantage of die large canister, however, is diat iu cloud 
reached the outer limits of the Uble much more quickly Uuin d« smaller 
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canister. The small canisler Ihus allowed longer event observations: in many 
cases reaching llie Tinal passive stage of tlie Iteavy gas event for which the 
cloud no lunger spreads primarily as a result of gravity forces. 

The obstacle array and canister size variations for the 18-iest series are 
presented in matrix form in Table I. Brief discussions of those tests are given 
below. For ease in the discussion, the cloud will be described as a flat, two- 
dimensional flgure with the "vertical" axis and direction indicated by the line 
running through the canister centre "up and down" on the photograph and the 
"horizontal" axis and direction taken as running "left and right" through the 
canister centre. Start tiiiKS for the releases shown in these figures arc given 
in Table 2 in hours, minutes, seconds and 1/30 second incrcttKnts corresponding 
to the lime signature shown in the individual photographs. 

Test Series 8 and 9 were run with no obstacles and with the small and 
large canisters lespeclivcly. In both cases, the heavy gas cloud spread in a 
circular manner as expected with the radius of the cloud increasing with time. 
The comparison of these clouds at an equal ume after release is shown in 
Figure 4. It should be noted that the clouds are nearly perfect circlci 
indicating that the experimental surface is very level and that there is no major 
influence by the release mechanism on the cloud dispersion. The upper 
photographs in Figure 4 are of the cloud M 0.33 seconds after release and the 
lower photographs are from 1.00 seconds after release. As expected, the cloud 
from the large canister has a larger radius than the cloud from the small 
canister after an equal time interval. It is also interesting to note the very 
visible outer "halo" at the cloud front and the thinner cloud density within die 
ring. (The bright inner circle is due to reflection firom the release mechanism 
pad.) 

Scries I, 2, S, 6 and 10 were run with full coverage (sheets of blocks both 
to the north and to the south of die canister) using both the small blocks and 
the medium blocks with die small and large canister releases. (Sec Table I for 
each setup for these series.) All of Uiese runs are characterized by a cloud 
pattern which initially begins as a circular shape but soon uansforms into a 
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diamond pailem wiUi apices ai the 'noilh-ioudi* ind 'cast-west" axes through 
the cylinder centre. An example of this flow pattern is given in Figure 5 for 
the large canister and medium blocks. The start ume for this scries was at 
03:40:00. As wiUi the "no blocks" cases, Uic large canister cloud moves more 
quickly Uian thai firom the small canister (photos of tmall canisler release not 
shown), 

Seriea 3, 4 and 7 were trials widi "half coverage", i.e. blocks only on the 
southern half of die study platform. (Sec Table 1 for conflguration details.) We 
sec dial the cloud patterns (Figure 6) generated by diis configuration arc half 
like die "no blocks* case and half like the "full coverage" case: a semicircle to 
die north and a half-diamond to die soudi. In the early life of Uic cloud, die 
spread is in a nearly perfect circular pattern widi cloud front movement 
through die blocks as rapid as dial in die block-fiee half. The most surprising 
result of dicse trials is dial the rate of spread of die cloud front is die tame 
along the circle radii as it is along die apex of dc diamond in die soudiem 
half. This similar rate of spread is found to occur regardless of die canister 
or block size. Tbe indicauon is diat, for diis half-coverage configuradon, a 
heavy gas cloud will flow along an unobstructed corridor between obstacles at 
the same rate as it would if no obstacles were prescnL Further discussion of 
diis phenomenon wUl be given in die next secuoo. 

Series II and 12 were trials using half and full coverage of medium-size 
blocks widi 1:1 spacing and die large canister release. The difference between 
diesc two series and Series 4, 3 and 10 is diat die alternate rows of blocks arc 
staggered by one block widdi. Widi diis anangcaient for die full coverage 
case, the open corridors along die north-soudi axis no longer exists. The 
arrangement does, however, sull have an open cast-west axis. The inidal cloud 
spreads as a circle (Figure 7) and dien transforms into a hexagon pattern. The 
half-coverage configuration results in an arrangement with a half-hexagonal 
pattern to the south and semi-circle to die north (Figure 8). The point of 
interest in these scries is diat die apices of die hexagon in die southern half 
(and in the northern half, in the case of full coverage) can be seen to occur in 
a diagonal direction from the release area, a path which is free of intersections 
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with any blocks. The movement along the north-south axis, however, is now 
impeded by the blocks and thus is slowed. The distance travelled from the 
canister centre to the centre points of the hexagon faces is essentially 
equidistant in all directions once the figure is established. The hexagon is 
elongated along the east-west axis as the cloud moves along the corridor free 
of impediment as it did in the regular array case. 

Scries IS and 16 consist of a heavy gas spill with the large and small 
canisters for full coverage of medium blocks at a 1:3 spacing. Again, the 
spread begins as a circle (Figure 9) and remains nearly circular although there 
is the hint of a dianmnd shape in the later times with the large canister. The 
edges of the cloud front are seen to be "scalloped" around the obstacles 
indicating the flow pattern around the blocks. 

Series 17 and 18 were run with the large canister and medium blocks 
with 2 cm-spacing (3.S:1 spacing) for full and half coverage respectively. 
Because of the narrowness of the corridors between blocks, some of the gas 
"sloshes" over the blocks upon release and forms a secondary cloud above the 
blocks. The pattern is further complicated by a rebound of the gas from the 
blocks. This process is dramatically seen in a series of scenes taken of the 
release region with a hand-held video camera (Figure 10). The cloud is seen lo 
break over the blocks, initially clearing the floor of the canister. The gas 
then returns to All the release area and a second wave is formed. This process 
is repeated with ever-diminishing strength. A side view of the process from a 
wider angle is shown in Figures 11 and 12. Very evident from this series is 
the secondary cloud spreading over the top of the blocks. This upper level 
cloud is evident firom the overhead photographs (Figure 13) as a brightening and 
obscuration of the block tops. The horizontal pattern initially spreads as a 
diamond and likely continues lo spread in that pattern although the gas 
becomes difflcutt to see in the corridors after some lime. The cloud moving 
over the blocks hints it a four-pointed star pattern which evolves toward i 
hexagon. 
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For the half-coverage case (Figure 14), the clear-area semi-circle is 

present widi a diamond shape evolving in the corridors between the blocks. The 

initial overflow of the blocks is evident as a rectangular bright area which 

'quickly dissipates. Since the rebounding material has open terrain ahead of it, 

there is no back-and-forth slosh as there was with the complete coverage. 

Series 13 and 14 were trials using simple solid barriers. In one, a pair of 
barriers were angled at 45 degrees and spaced so that a gap existed between 
them. In the other, a single barrier was placed horizontally. The cloud front 
was unaffected in all directions except where it touched the barrier which 
stopped the cloud front. The cloud wrapped around the barriers and would 
eventually fill in the space behind. Figure 15 shows these releases about a 
second after release. 

4. ANALYSIS OF FLOW VISUALIZATION RESULTS 

The flow visualization experiments provided a good record for analyzing 
the movement of the heavy gas cloud after iu release. By reducing the 
videoupc record to a series of photographs, the horizontal dimensions of die 
cloud could be determined and a time/distance history of die cloud front under 
a variety of release conditions ascertained. This was accomplished by measuring 
the distance on an enlargement of the photograph of various poinu along die 
cloud front from the canister centre and scaling the result by the known 
distance between die horizontal tapes above and below the canister. 

The error in die measurement of points on die cloud front was ±0.5 nun 
from the photographs. The scale factor, between the photograph and die 
model, ranged from 22 to 2S. Thus, die actual error in the placement of die 
cloud front is 1.1 cm to 1.3 cm. The clock used resolved time into increments 
of 1/30 of a second. 

The distance of the cloud front from the canister centre as a function of 
elapsed time is shown in Figure 16 for both the large and small canisters 
under "no blocks" configurations. As expected,, flow from die large canister 
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spread at a faster rate than from the small canister. The solid lines represent 
Ihe current MOE model with coefricieni equal to unity. Modelling, including 
the dennitlon of Ihe coefricieni, is discussed in detail in the next section. 

In this paper, only the 1:1 regular block spacing has been analyzed for 
cloud spread. For the large canister, the cloud front distance as a function of 
elapsed time is shown in Figures 17 and 18 for a variety of block size and 
coverage conflgurations. Figure 17 shows the distance measured along Ihe axes 
(90 degrees) and Figure 18 shows die distance as measured at a 4S degree angle 
through die blocks. Similarly, Figure 19 (90 degrees) and Figure 20 (43 
degrees) show measund distances for Uie small canister under several block 
configurations. 

When compared with the configuralions with blocks, die frontal distances 
for large canister releases along the axes for the "no blocks" connguration are 
less than distances measured with blocks present as full coverage (Figure 17). 
In other words, the cloud front moves faster with blocks present Uian widioul, 
at least along die corridors, and this holds for full coverage and no-blocki 
conflgurations only. Despite some scatter In die dau, diere is a distinct 
pattern of die distances aligning in relation to the block conflguration. The 
distance travelled after a given time, i.e. die cloud spread, is seen to increase 
in the following sequence of block patterns: no blocks, medium blocks-half 
coverage, medium blocks-full coverage, small blocks-full coverage. The spread 
of Ihe cloud from die small canister, however, shows only t slight 
differentiaUon between conflgurations with die 'no blocks" spread being 
marginally faster (Figure 19). The velocity difference between conTigurauons 
may be due to the different combinations in the source size, obstacle size and 
obstacle spacing but how die process evolves has not yet been fully deteimined. 
ll is interesting lo note diat, in the half-coverage cases, die rate of spread 
along the north-south axes, through blocks on soudiem half and open in 
nonhcrn half, is virtually identical. For full block and no block coverage, the 
reason is obvious, Ihe symmetry of obsucle placement. However, why Uiis 
holds also for the half coverage cases is also not yet understood. 
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The distances travelled along die 43 degree padis (Figures 18 and 20) do 
not show a distinct pattern among block cdnHgurations, during the initial suges 
of die spill. In later stages, however, die rate of spread was higher for die 
"no-blocks" conflguration. This difference in spreading rate seems to 
distinguish a transition from a circular spreading pattern (initial suge) to a 
diamond spreading pattern (next suge). The transition occurs earlier for a 
large canister release (0.6 - 1,0 seconds) dian for die small canister release 
(1.0 - 1.3 secciids). This was expected because of die higher presaure head 
afforded by die large canister. 

5. MODELLING OF HORIZONTAL CLOUD SPREAD 

The atmospheric dispersion of heavy gases cannot be described fully by 
general gas dispersion mediods. Madiematical models which treat die gravity 
spreading of a heavy gas cloud and die entrainmenl of ambient air into die 
cloud have been developed by a number of investigators over die past two 
decades. An excellent review was presented by Havens (1982). The modelling of 
an instantaneous heavy gas release is generally considered to be comprised of 
ditee phases. The first phase is die "blow-up" phase where die gas is initially 
released into the environmenL The second phase is the "initial spreading" 
phase where the cloud moves as a slumping column of nearly uniform gas which 
mixes via turbulent entrainment of ambient air. The second phase ends and Ihe 
third phase begins when the density difference between the cloud and ambient 
becomes sufficiendy small. In diis third phase, die gas cloud is sufflciendy 
dilute dial iu dispersion may be treated by die traditional puff dispersion 
dieoties. 

We consider here die initial suges of die second phase. It is assumed 
dial die result of die gas release is a cylindrical cloud of known dimensions. 
Using a simple box model approach to heavy gas dispersion (MOB, 1983), die 
gas slumping phase is assumed to be controlled by graviutional effects. 
Assuming a constant cloud volume during diis phase gives the following 
relationship for the rate of spread of die cloud under still air conditions: 
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(3) 



where r is Ihe cloud radius at time I. r^ is llie initial cloud radius at 1=0, C is 
a constant, A6 = 8 - 6^, 5 and 8^ are the gas and air densities respectively, 
Vq is the initial gas volume, and g is Ihe gravitational constant. 

In theory, Ihe constant C could be either 1 or 1.414 (Hoult, 1972), but 
experimental results have shown C=l to be the belter choice (van Ulden, 1974; 
de Nevera. 1984). We next consider the value of C in light of our data. 

Applying equation (3), with C=l, to the experimental data from the "no 
blocks" case (presented in Figure 16) results in a slight under prediction of 
the distance to the cloud front. (Note the lines in Figure 16 refer to the 
model values). The model consistently underpredicts the cloud radius for both 
canisters. In order to improve the model fit, the constant C was determined 
empirically through (3) for each time segment and the resulting values averaged. 
Under Ihe "no block" conriguratioo for the large canister, the average value of 
C was found to be 1.19. When this value for C is used in (3), the lower curve 
on Figures 17 results. The upper set of points in Figure 17 are the resulu of 
the "small block-full coverage* scenario. The value for C calculated for this 
configuration is 1.75, and a plot of (3) using this value of C is shown as the 
upper solid curve in Figure 17. A value of C-1.34 which ii an avenge of all 
values of C from all large canister releases with blocks gives a line (noc 
shown) midway between the two extremes. 

For the small canister releases, the value of C determined from several 
different configurations ranges between 1.24 and 1.27. Only the plot of (3) 
with C>1.27 is shown on Figure 19. 

From these physical model experiments, the "constant" C in equation (3) is 
found to vary with release size and, for Ihe large release, with the obstacle 
conflguralion. We next consider whether Uiis constant can be related to block 
size or spacing. 
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In diose releases wiUi a regular 1:1 block array, the heavy gas cloud was 
observed to begin as a circular cloud and then take on a diamond shape as it 
moved through die blocks. This flow regime shows the portion of the cloud 
moving down the unobstructed corridors between the blocks along the two 
horizontal axes to be spreading die quickest while the flow on the 4S° 
padiway moves slowest after an initial period of equal speed. From Ihe 
geometry of Ihe final configuniion, a simple model of the flow regime is 
proposed to be: 



rj - r f(t) 



(4) 



where r^ is the distance of spread from the canister centre, r is the distance 
calculated by (3) down die axes, f(t) is the shape factor u * fiinction of 

dme L Note f(l) is defined as: 

f(t) - lain 45/sin(135-a)l[e-*(sin(135-o)/sin45) + (l-e"*)! (5) 

where a is die angle (in degrees) of die radius from die horizontal axis and S 
is a factor which is equal to 1 for blocks and for no blocks. Note dial for B 
= 0, f(l) -» 1, and equation (4) becomes r^ - r, i.e. Ihe equation of a 
circular cloud front Widi B = I, as t -> — , equation (5) reduces to 
sin(4S)/sin(13S-a). Equation (4) dien describes a linear cloud front connecting 
the two maximum radii at o^^" and 90°. Figures 18 and 20 show die results of 
dii» model for die 4S° radius for die large and small canisters respectively. 
In Figure 18, die distance wu calculated using 01.S4, die value which gave 
the best fit for all die block thita. 

The mathematical model for heavy gas cloud radius (gas cloud spread) as a 
function of time given by (4) provides a good fit to die measured cloud 
dimensions under calm winds for dispersion diough a regular array of blocks as 
well as die "no blocks" situation. This model is similar in form to the 
graviutional spread models proposed by many audiors (e.g. van Ulden, 1974; 
MOE, 1983; de Nevers, 1984) widi die "constant" term, generally taken as I, 
increased to between 1.19 and 1.7S depending upon die obstacle configuration. 
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This very specific approach of ruling data niusi be generalized to incoqwrale 
other lypes of complex leirain. As more dala become available, Iwo approaches 
may be considered. Firstly, dala from several dirferenl terrain types may be 
fllled lo a speciHc nradel, where for each terrain type a different set of 
empirical constants are recotnroended. Tliis leaves the end user to make a 
judgement on the terrain type. A second approach is to determine a set of 
empirical constants averaged over all terrains tested, and including the error 
battd. The end user can than use the error bands in making judgements on the 
cloud path. In subsequent work we shall look at both approaches in light of 
wind tunnel data on heavy gas conceniiation mcasureiiKnts. 

6. CONCLUSIONS AND FUTURE WORK 

The flow visualizauon experiments for an instantaneous heavy gas released 
in an array of regularly spaced obstacles have resulted in a number of 
interesting spread patterns. The analysis of the visual record of these releases 
has given a cloud spread history for each of the various configurations wliich 
has been used to test a basic mathematical model of cloud radius under 
gravilalional spreading. The results of thai comparison indicate that die models 
nt the observations well when an empirical constant is introduced. 

The next step in this study is to measure the gaseous concentration at a 
number of locations and lo expand the gas spread model lo predict 
concentrations. In addidon, other release conditions will be itudied including 
the non-zero wind case.. 
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TABLE I 
TEST SERIES CONFIGURATIONS 



Test Series Canisler 
Number Size 



Block 
Size 



Block Block 
Spacing Coverage 



I 

2 

3 

4 

5 

6 

7 

I 

9 

10 

II 

12 

13 

f4 

13 

IS 

i7 
It 



S 
L 
L 



S 

L 
L 
L 
L 
L 
L 
L 
S 
L 
L 



S 
S 

s 

M 
M 
M 
M 



M 
M 
M 



U 
M 

M 
M 



1:1 
1:1 
1:1 
1:1 
1:1 
1:1 
1:1 



1:1 
1:1 
1:1 



1:3 

1:3 

3.75:1 

3.75:1 



Full 

FiiU 

Italf 

Half 

Full 

Full 

Half 

None 

None 

Full 

Full Staggered 

Half Staggered 

2 Barriers, 45° 

1 Barrier 

Full 

Full 

Full 

Full 



TABLE 2 
START TIME FOR EXAMPLE RELEASES IN PHOTO SEQUENCES 



Figure Number 



2 
4* 

4k 

i 

6 
7 
I 

9 

10 

11 

12 

13 

14 

ISt 

I3b 



Start Time 



01:23:36:05 
00:03:24:05 
00:02:42:20 
00:03:40:00 
00:03:57:15 
00:08:37:20 
00:09:17:05 
00:08:18:00 

00:04:51:11 
00:04:51:11 
00:04:14:20 
00:05:44:15 



Legend: Canister Size: S - small; L - large 
Block Size: S - small; M - medium 
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nCURE 3 
UNIFORM TERRAIN IN THE WIND TUNNEL 
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RGURE 5 

OVERHEAD VIEW OF HEAVY GAS RELEASE: LARGE CANISTER, 

MEDIUM BLOCKS, 1:1 SPACING, FULL COVERAGE 
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RGURE 6 

OVERHEAD VIEW OF HEAVY GAS RELEASE: LARGE CANISTER. 
MEDIUM BLOCKS, 1:1 SPACING, HALF COVERAGE 
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RGURE 7 

OVERHEAD VIEW OF HEAVY GAS RELEASE: LARGE CANISTER. 
MEDIUM BLOCKS, 1:1 STAGGERED SPACING. FULL COVERAGE 
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RGURE 8 
OVERHEAD VIEW OF HEAVY GAS RELEASE: LARGE CANISTER. 
MEDIUM BLOCKS. 1:1 STAGGERED SPACING. HALF COVERAGE 
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RGURE 9 

OVERHEAD VIEW OF HEAVY GAS RELEASE: LARGE CANISTER. 
MEDIUM BLOCKS, 1:3 SPACING, FULL COVERAGE 
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RGURE 10 

RELEASE AREA FOR HEAVY GAS RELEASE: LARGE CANISTER. 
MEDIUM BLOCKS, 3.5:1 SPACING, FULL COVERAGE 
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RGURE 11 

SIDE VIEW OF HEAVY GAS RELEASE - EARLY STAGES: LARGE CANISTER. 

MEDIUM BLOCKS. 3.5:1 SPACING, FULL COVERAGE 




FIGURE 12 

SIDE VIEW OF HEAVY GAS RELEASE - LATE STAGES: LARGE CANISTER. 
MEDIUM BLOCKS. 3.5:1 SPACING, FULL COVERAGE 
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RGURE 13 

OVERHEAD VIEW OF HEAVY GAS RELEASE: LARGE CANISTER. 
MEDIUM BLOCKS. 3.5:1 SPACING. FULL COVERAGE 




RGURE 14 

OVERHEAD VIEW OF HEAVY GAS RELEASE: LARGE CANISTER. 
MEDIUM BLOCKS. 3.5:1 SPACING. HALF COVERAGE 





RGURE 15 

OVERHEAD VIEW OF HEAVY GAS RELEASE: LARGE CANISTER 
A SOUD BARRIER, B TWO BARRIERS WITW GAP 
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FIGURE 16 
DISTANCE vs TIME: NO BLOCKS 
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FIGURE 17 
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DISTANCE vs TIME: LARGE CANNISTER 90deg 
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FIGURE IB 
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DISTANCE vs TIME: LARGE CANNISTER 45deg 
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FIGURE 19 
DISTANCE vs TIME: SMALL CANNISTER 90deg 
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RGURE 20 
DISTANCE vs TIME: SMALL CANNISTER 45deg 



